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Abstract. We explore the viability of the unification of 
BL Lacs and FR I radio galaxies by comparing the core 
emission of radio galaxies with those of BL Lacs of simi- 
lar extended radio power, taking advantage of the newly 
measured optical nuclear luminosity of FR I sources. The 
spectral properties of complete samples are also studied 
in the radio-optical luminosity plane: starting from the 
Spectral Energy Distribution (SED) of BL Lacs, we cal- 
culate the predicted luminosity of FR I nuclei in the frame 
of a simple one-zone model, by properly taking into ac- 
count the relativistic transformations. We find that the 
bulk Lorentz factors required by the spread in the ob- 
served luminosities in all bands are significantly smaller 
than those implied by other, both observational and the- 
oretical, considerations. This discrepancy is also reflected 
in the fact that FR I nuclei are over-luminous by a fac- 
tor of 10-10 4 , with respect to the predictions, both in the 
radio and in the optical band. 

In order to reconcile these results with the unification 
scheme, velocity structures in the jet are suggested, where 
a fast spine is surrounded by a slow (but still relativistic) 
layer so that the emission at different angles is dominated 
by different velocity components: the fast one dominates 
the emission in BL Lacs while the slow layer dominates 
the emission in misaligned objects. Furthermore for the 
lowest luminosity BL Lacs it has to be also postulated that 
their beaming factor in the radio band is lower than in the 
optical (and X-ray), as would result from deceleration of 
the jet. 

The self-consistency of the unification model there- 
fore requires that both intrinsic differences in the SED 
and different beaming properties play a substantial role 
in characterizing the phenomenology of these sources. 

Key words: Galaxies: active; Galaxies: BL Lac objects: 
general; Galaxies: jets; Galaxies: nuclei 



1. Introduction 

Unification models adduce the main differences between 
the observed properties of different classes of AGNs to 
the anisotropy of the radiation emitted by the active nu- 
cleus (see Antonucci 1992 and Urry & Padovani 1995| for 
reviews). In particular, for low luminosity radio- loud ob- 
jects, namely BL Lacs and FR I radio galaxies (Fanaroff 
& Riley 1974), it is believed that this effect is mainly due 
to relativistic beaming. In fact, there is growing evidence 
that obscuration does not play a significant role in these 
objects, contrary to other classes of AGNs. This is in- 
dicated by optical (Chiaberge et al. 1999 hereafter Pa- 
per I), radio (Henkel et al. 1998), and X-ray information 
(e.g. Fabbiano et al 1984| Worral and Birkinshaw 1994. 
Trussoni et al. 1999 ). Within this scenario, the emission 
from the inner regions of a relativistic jet dominates the 
observed radiation in BL Lacs, while in FR I, whose jet is 
observed at larger angles with respect to the line of sight, 
this component is strongly debeamed. Evidence for this 
unification scheme includes the power and morphology of 
the exten ded r adio emission of BL La cs (e.g. Anton ucci fc 
Ulvestad |1985| , Kollgaard et al. |l992| , Murphy et al. |l993|) 
and the prop erties of their hos t gala xies (e.g. Ulrich 1989 , 
Stickel et al. 1991 , Urry et al. 1999 ), which are similar to 
those of FR I. Furthermore, there is a quantitative agree- 
ment among the amount of beaming required by different 
observational properties (e.g. Ghisellini et al. 1993), the 
number densities and luminosity functions of the parent 
and beamed populations in different bands (e.g. Urry & 
Padovani |1995| , Celotti et al. |1993| ) and the comparison of 
the radio core emission of beamed and unbeamed objects 
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with similar total radio power (Kollgaard et al. 1996). 

Despite this global agreement, it should be stressed 
that beaming factors inferred from the broad band spec- 
tral properties of blazars, more specifically super luminal 
motions, transparency to the 7-ray emission, shape of the 
SED and time-lags among variations at different frequen- 
cies, are significantly and systematically larger than those 
sugge sted by radio l umino sity data (Dondi fc G hisellini 
1995j Ghisellini et al. |l998j Tavecchio et al. |l998|) . 
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Thanks to the Hubble Space Telescope (HST), faint 
optical nuclear components have been recently detected 



in FR I galaxies (Chiaberge et al. 1999). A strong linear 
correlation is found between this optical and the radio core 
emission which strongly argues for a common non-thermal 
origin. This suggests that the optical cores can be identi- 
fied with synchrotron radiation produced in a relativistic 
jet, qualitatively supporting the unifying model for FR I 
and BL Lacs. 

These information offer a new possibility of verifying 
the unification scheme, by directly comparing the proper- 
ties of the optical and radio cores of radio galaxies with 
their putative aligned (beamed) counterparts, analogously 
to the procedure followed for the radio cores. X-ray obser- 
vations also provide useful constraints to the nuclear emis- 



2.2. Radio and X-ray selected BL Lacs samples 

We consider both the complete sample of 34 radio selected 



sion of FR I sources (e.g. Hardcastle & Worrall 1999). 

The main advantage of using multifrcquency data is 
the possibility of directly comparing the full broad band 
spectral distributions of these two classes of sources and 
eventually shed light on the apparent discrepancy in the 
Lorentz factors inferred from different approaches. 

The paper is organized as follows. The (complete) sam- 
ples of BL Lacs and radio galaxies are presented in Sect. 0. 
In Sect. H we compare separately the core radio and opti- 
cal emission of beamed and unbeamed objects with similar 
extended radio power. From this we infer the Lorentz fac- 
tors requested by the unification scheme within the sim- 
plest scenario in which the radiation is emitted by a single 
uniform region of the relativistic jet. In Sect. |i| the ra- 
dio and optical data are considered together and, starting 
from the observed SED of BL Lacs, we derive the expected 
properties of the nuclear emission of FR I, by taking into 
account the spectral dependence of the relativistic trans- 
formations. As the single-region picture does not account 
for the observed properties, in Sect. ||we explore a (sim- 
ple) alternative scenario and test it also against the X-ray 
information. Summary and conclusions are presented in 
Sect, g 

2. The samples 

2. 1 . FR I radio galaxies 

Our complete sample of radio galaxies comprises all the 
FR I sourc es belonging to the 3CR catalogue (Spinrad 
et al. 1985), morphologically identified as FR I. The red- 
shifts of these objects span the range z = 0.0037 - 0.29, 
with a median value of z — 0.03, and the total radio lumi- 
nosities at 1.4 GHz are between 10 30 2 and 10 34 ' 2 erg s _1 
Hz -1 (H = 75 km s -1 Mpc -1 and go = 0.5 are adopted 
hereafter). We exclude from the original sample the pe- 
culiar object 3C 386, as discussed in Paper I. The optical 
and radio data are from Paper I, while the X-ray ones 



BL Lacs derived from the Uy catalog (Stickel et al. 1991 



Kuhr et al. 1981), and the BL Lac sample selected from 
the Einstein Slew survey (Elvis et al. 1992, Perlman et al. 
1996), which comprises 48 objects, and it is nearly com- 



plete. The extended radio power (at 1.4 GHz) L ext [erg 
s" 1 Hz" 1 ] span s the ranges 10 30 1 - 10 33 8 (Uy BL Lacs, 
Kollgaard et al. |1996| ) and 10 29 1 - 10 33 4 ( Slew s urvey BL 
Lacs, Kollgaard et al. 1996, Perlman et al. [1996 ) ; the red- 
shifts are between 0.049 and 1.048 (median z = 0.501) and 
between 0.031 and 0.513 (median z — 0.188) for the two 
samples, respectively (the redshifts of the Slew B L Lac s 
are taken from the data collected by Fossati et al. 1998 ). 



Instead of classifying BL Lacs according to their selec- 
tion spectral band, in the following we adopt the defini- 
tions of high and low energy peaked BL Lacs (HBL and 
LBL respectively) , which are based on the position of the 
(synchrotron) emission peak in the spectrum and there- 
fore more indicative of the physical characteristics of the 
objects (Giommi & Padovani 1994, Fossati et al. 1998). Of 
the 34 objects belonging to the 1 Jy sample, 32 are classi- 
fied as LBL and 2 as HBL, while of the 48 X-ray selected 
BL Lacs, 40 are HBL and 8 are LBL. 

The spectral data fo r bot h samples of BL Lacs are 
taken from Fossati et al. 1998 



from Hardcastle & Worral (1999). The optical data are 
extrapolated to the V band using a spectral index a a = 1 
{F v (xu~ a ). 



3. Core versus extended luminosity 

According to the unification models the beamed and un- 
beamed populations must cover the same range of ex- 
tended luminosity, as this is considered to be isotropic. On 
the contrary, emission from the core is affected by beam- 
ing: radio galaxies should have a fainter central compo- 
nent, whose intensity would depend on the Doppler factor 
5 = [r(l - f3cos9)}-\ where T = (1 - /3 2 )~ 1/2 , (3c is the 
bulk velocity of the emitting plasma and 9 the angle be- 
tween the direction of the jet and the line of sight. The 
transformation law for the specific flux density is in fact 
F v — 8 p+a Fl, where the primed quantity refers to the co- 
moving frame, a is the local spectral index, p = 2 for a 
continuous jet and p = 3 for a moving sphere. 

Therefore the comparison of the core emission of 
beamed objects and their parent population with simi- 
lar extended emission provides a direct estimate of the 
Lorentz factor of the radiating plasma, if the typical ob- 
serving angles are known. With this aim and similarly to 
what has been done in the radio band (e.g. Kollgaard et 
al. 1996), we plot the optical V band luminosity (L ) vs 
the extended radio luminosity at 1.4 GHz (L ext ) for the 
three samples (Fig. Q). 

First we should note that the HBL objects do not fully 
share the range of extended radio power of the 3CR ra- 
dio galaxies (the HBL total luminosities are in fact more 
similar to the objects belonging to the B2 sample of low 
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Fig. 1. Optical core luminosity (V band) versus radio ex- 
tended luminosity at 1.4 GHz for FR I (circles), HBLs 
(squares) and LBLs (triangles). The grey scale refers to 
the three bins of extended radio power. The dashed line 
represents the linear fit to the FR I sample, having ex- 
cluded the most aligned sources, here marked with crosses 
(see text). The range of extended power covered by the B2 
sample of radio galaxies is also indicated. 



power radio galaxies). Conversely, L ext of LBL well match 
those FR I of the 3CR catalog. 

Also notice that the regions occupied by the two sam- 
ples of BL Lacs appear to be continuously connected, the 
lower radio power BL Lacs (which are HBLs) and the 
higher radio power ones (LBLs) having an optical luminos- 
ity which weakly increases for increasing extended lumi- 
nosity. Because of this trend, in order to compare sources 
with the same L ext we have sub-divided the samples into 
three bins, namely: \ogL ext [erg s _1 Hz -1 ] < 31.5, \ogL ext 
between 31.5 and 32.5, and \ogL ext > 32.5. 

We thus calculate the median values of the observed 
nuclear luminosity of FR I and BL Lacs in each interval 
of extended power. BL Lacs are on average 4 orders of 
magnitude brighter than FR I cores. We can assume that 
BL Lacs are observed [] at ~ 1/r and FR Is at = 60°: 
in fact, for an isotropic distribution of objects, = 60° 
corresponds to the median angle if, as it is in the case of 
FR I, the scatter in the optical luminosity is dominated 
by relativistic beaming. Bulk Lorentz factors r ~ 4 for 
the case of an emitting sphere and ~ 6 for a continuous 
jet are required in order to account for the different core 
luminosities of FR I and BL Lacs in each bin of extended 
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Fig. 2. The regions occupied by the three samples in the 
optical luminosity versus extended radio luminosity plane, 
as for Fig. |l|. The dashed lines indicate the correlation 
found between these two quantities when shifted by beam- 
ing effects for the values of the bulk Lorentz factor marked 
on the left. 



power. An optical spectral index a = 1 is assumed for all 
sources (independent of beaming) . 

An alternative method to estimate T relies on the fact 
that, for a randomly oriented sample, the best fit regres- 
sion line of a luminosity distribution corresponds to the 
behavior of sources observed at ~ 60° , once the most core 
dominated objects are excluded (Kollgaard et al. 



1996|) 



At the particular angle ~ 1/r, one obtains 5 = T. 



We thus determine the best fit regression of FR I in the 
L ext — L a plane, after excluding from the sample the 5 
objects in which optical jets are detected. These sources, 
in fact, have the most luminous optical cores, are among 
the most core dominated objects in the radio band, and 
their radio jets are shorter, indicating that they are point- 
ing towards the observer (Sparks et al. 1995). Interest- 
ingly, we obtain that there is a remarkable correlation 
(P > 99.9%) between \ogL ext and logL Q , among the re- 
maining 20 "highly misoriented" objects, although with a 
slope (~ 0.6) marginally steeper than the correlation be- 
tween L ext and core radio luminosity (L r ) found by Gio- 
vannini et al. ( 198§D for a larger sample of radio galaxies. 
In Fig. ^| we show the regions in which the three samples 
are located in the L ext — L Q plane, and the dashed lines 
represent the "beamed" FR I population as observed un- 
der an angle ~ 1/r in the case of p = 3. Also with this 
method r ~ 5 (~ 7) are required to displace the FR I 
to the regions occupied by both HBL and LBL for p = 3 
(p = 2). 
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Fig. 3. Radio core luminosity (at 5 GHz) versus radio ex- 
tended luminosity at 1.4 GHz for FR I (circles), X-ray 
selected (squares) and radio-selected BL Lacs (triangles). 
The dashed line is the correlation between these two quan- 
tities found for a larger sample of galaxies by Giovannini 
et al. ( 1988 ) and converted to 1.4 GHz using a r = 0.7. 



Let us now consider L r (at 5 GHz) versus L ext (Fig. |J), 
analogously to what is shown by Kollgaard et al. ( 1996 ) for 



a larger sample of radio galaxies (which also includes our 
objects). The typical radio core luminosities of HBL and 
LBL are significantly different, the latter objects being on 
average about one order of magnitude more luminous than 
the former ones. Conversely, as we have already pointed 
out, no substantial difference between the two classes is 
found in the case of L . 

These results have been initially attributed to a dif- 
ferent amount of beaming for X-ray and radio-selected 
BL Lacs (i.e. different angle of sight and/or different jet 
velocities^) while more recently a consistent picture has 
emerged where this diversity can be accounted for by the 
different shape of their intrinsic SED (e.g. Padovani 

Giommi & Padovani 



1992 



Ghisellini & Maraschi 1989 



Fossati et al. 1998Q . The role of these two scenarios will 



be further explored in the next section, through the com- 
parison of the SED of both types of BL Lacs with their 
parents. 

We conclude that the Lorentz factors inferred from the 
comparison of the radio, but also optical emission of FR I 
and BL Lacs, are consistent with those previously esti- 



2 If so the inferred Lorentz factors, relative to our sample, 
are T = 4 (3) for HBL Y = 10 (5) for LBL, in the case of p=2 
(P=3) 
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Fig. 4. BL Lacs and FR I radio galaxies in the L r — L a 
plane. Empty symbols are objects with no data on their 
extended radio power, filled symbols and grey scale are as 
in Fig. [l]. The dashed line is the radio-optical core corre- 
lation (Paper I). The range of core luminosity of the B2 
radio galaxies is also reported. Notice that the B2 cores 
could at most extend this correlation by ~ 1 order of mag- 
nitude towards lower luminosities. 



mated from the statistics of these sources within the uni- 
fying scheme. However, as already mentioned, such val- 
ues are significantly and systematically lower than those 
required by other independent means, such as superlumi- 
nal motions and high energy spectral constraints (fit to 
the overall SED and ti me-lags) in bo th LB Ls and HBLs. 
(Mar aschi et al. |l992j, Sikora et al. |1994 Celotti et al. 



1998, Tavecchio et al. 1998). These latter methods require 



a value of the Doppler factor S in the range 15-20 for the 
region emitting most of the radiation in both HBLs and 
LBLs. The need for high degrees of beaming will consti- 
tute a crucial point in the following. 



1994, 4. FR I and BL Lac in the L a and L r plane 



Since for the first time multifrequency data are available 
also for the nucleus of radio galaxies, we can now directly 
compare the spectral properties of beamed objects and 
their parent population: this new approach can thus com- 
bine information from the SED of BL Lacs with those in- 
ferred from their relation with FR I. In particular in Paper 
I and in Chiaberge et al. 1999 we showed how the loca- 



tion of sources in the L D and L r plane represents a very 
useful tool to discuss their nuclear properties. In Fig. |J 
we show the optical vs radio core luminosity for the three 
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samples. The dashed line represents the (almost linear) 
correlation found between these two quantities among the 
FR I sources (Paper I). Radio galaxies, HBL and LBL 
occupy different regions of this plane: LBLs are located 
only marginally above the continuation of this correla- 
tion, while HBL are ~ 2 order of magnitude brighter in 
the optical with respect to other objects for a given radio 
luminosity. 

In order to determine how beaming affects the ob- 
served luminosities and thus how objects could be con- 
nected in this plane, we consider the SED of BL Lacs, ob- 
servationally much better determined, and calculate the 
observed spectrum of the misaligned objects, by taking 
into account relativistic transformations. 

In fact, an important point, previously neglected, is 
that these transformations depend on the spectral index 
in the band considered, which in itself might change as 
a function of the degree of beaming. Therefore, in order 
to correctly de-beam the SED of BL Lacs, a continuous 
representation of it and an estimate of the bulk Lorentz 
factor of the emitting region are needed (it is again as- 
sumed 9 = 1/r). While any continuous description of 
the SED and typical Lorentz factors can be used, we de- 
rive both of them by adopting a homogeneous synchrotron 
self-Compton emission model to reproduce the observed 
SED s (e.g. Chiaberge & G hiscll ini |1999| , Ghisellini et al. 
199S| , Mastichiadis & Kirk |l997j ). This approach has the 



advantage of considering both the emission and dynamical 
(F) properties self-consistently. The bulk velocities ob- 
tained in this way are fully compatible with those inferred 
from the already mentioned constraints (Sect. |^). 

Let us firstly consider single objects for which the SED 
is well sampled, namely Mkn 421 (a typical HBL) and PKS 
0735+178 (a typical LBL). We model the observed SED 
as explained, derive the value of T(— S) for the two sources 
and then calculate their corresponding observed SEDs for 
different orientations. Clearly the net effect of debeaming 
is a "shift" of the SED towards lower luminosities and 
energies (see Fig. ||). 

Notice that, as the model is appropriate for the op- 
tically thin part of the spectrum, in order to account 
for the radio emission, which necessarily has to be pro- 
duced on larger scales, we linearly extrapolate the fit from 
the infrared-mm spectral region. However, at an angle of 
9 = 60° and for the Lorentz factors derived from the 
model, r ~ 15 — 20, the observed (debeamed) radiation at 
5 GHz corresponds to what is seen in BL Lacs at far in- 
frared frequencies (respectively ~ 500 — 300jum, see Fig. ||) 
and therefore the debeamed points in Fig. o represent the 
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Fig. 5. Spectral energy distributions of Mkn 421 (lower 
panel) and PKS 0735+178 (top panel) and debeamed 
SED for different viewing angles in the case of a sin- 
gle emitting component. The Lorentz factors inferred for 
the two sources are r — 20 and T — 16, respectively. 
For comparison, we report (empty circles) the radio, IR 
and optical (HS T) an d X-ray (ROSAT) data for 3C 264 
(Capetti et al. 2000). The (non-simultaneous) data for 
PKS 0735+178 are taken from the literature (NED). The 
(quasi-simultaneous) data for Mkn 421 are from Macomb 
et al. (|l99l). 



These models, in fact, satisfy time variability constraints, 
assume continuous injection of particles, radiative and adia- 
batic cooling, 7-7 collisions and pair production. The contin- 
uous curves shown in Fig. 5 are then not only interpolating 
curves, but physically possible fits to the data. 



correct predicted luminosities of the BL Lac component 
at 5 GHz. 

The resulting debeamed optical and radio luminosities 
are reproduced in Fig. |[ The dash-dotted lines represent 
"debeaming trails" and the filled circles the calculated de- 
beamed luminosities for 9 — 1/T (i.e. the BL Lac itself), 
10°, 30° and 60°. Most noticeably, for 9 = 60° - which is 
the mean angle of sight for the misaligned population - 
the BL Lac component is about four orders of magnitude 
below the radio galaxy region in the optical, and two/four 
in the radio band. 

While equivalently incompatible with the FR I popu- 
lation, the HBL and LBL move on different trails. This is 
due to the different shape of their SED (see Fig. |J), and 
in particular to the position of the synchrotron peak fre- 
quency: if - for increasing values of 9 - in the rest frame 
the peak overcomes the optical band, the spectral index 
steepens and the optical flux drops more rapidly than the 
radio one. 

Another remarkable result is that the debeaming trail 
of the HBL does not even cross the region occupied by 
radio galaxies in the L r — L Q plane. As this might be a se- 
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Fig. 6. Debeaming trails in the radio-optical luminosity 
plane for Mkn 421 and PKS 0735+178, in the frame of a 
single emitting region model. The filled circles correspond 
to the predicted luminosities of objects at different viewing 
angles. Top to bottom: 6 = 1/T, 10°, 30° and 60°. 
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Fig. 7. Debeaming trails in the radio-optical luminosity 
plane for average BL Lacs SED, in the frame of a single 
emitting component model. The filled circles correspond 
to the predicted optical-radio luminosity for different an- 
gles of sight (top to bottom: 9 = 1/T, 10°, 30° and 60°). 



rious problem for the unified scheme, we further examine 
this issue. In particular we closely examine the effect of 
the spectral shape and its relation with the intrinsic lumi- 
nosity by considering three different SED, which represent 
the whole family of BL Lacs, from HBLs to LBLs. Qln Fig. 
we plot the resulting trails: once again, as in the cases of 
Mkn 421 and PKS 0735+178, the expected nuclear lumi- 
nosity is 10-10 4 times fainter than what observed in FR I 
and the debeaming trail for the lower luminosity object (a 
typical HBL) does not even cross the FR I region. Note 
that if the luminosity is indeed related to the shape of the 
SED, this discrepancy would exacerbate for even fainter 
BL Lacs. In fact for HBL the radio and optical spectral 
slopes can be considered constant as the viewing angle in- 
creases, resulting in a linear (one to one) debeaming trail, 
parallel to the FR I correlation. 

Summarizing: the radio and optical luminosities of BL 
Lacs and FR I are not consistent with the simplest pre- 
dictions of the unifying scheme, if a single emitting region 
is responsible for the different broad band spectral prop- 
erties of the beamed and parent populations. More specif- 
ically: i) Lorentz factors ~ 15, as derived from the high 
energy spectral properties, underestimate the predicted 



4 The three SED correspond to different bins of radio lumi- 
nosity (at 5 GHz) - which appears correlated with the bolomet- 
ric lu minos ity and the posit ion o f the peak frequency (Fossati 
et al. [l998|, Ghisellini et al. |l998|). 



emission from the parent population; ii) the relative ratio 
of radio to optical luminosity of HBL is inconsistent with 
the observed FR I spectra. In the next section we discuss 
and test a possible solution to these discrepancies. 

5. A jet velocity structure 

We have shown that the high bulk Lorentz factors required 
by the emission models of BL Lacs imply that if such ob- 
jects are observed at 60° the resulting spectral properties 
are not compatible with what is observed in the nuclei of 
radio galaxies. And indeed the previous comparison of the 
core emission of FR I and BL Lacs (see Sect. ||) led to lower 
values of T. How can these results be reconciled within the 
unifying scenario? A possible and plausible effect, which 
could account for this discrepancy, is provided by the ex- 
istence of a distribution in the bulk velocity of the flow, 
with the emission from plasma moving at different speeds 
dominating the flux observed at different viewing angles. 

Let us consider this hypothesis in the frame of the uni- 
fication scheme and examine the simplest case, i.e. a model 
with two axisymmetric components having the same in- 
trinsic luminosity and spectra. In other words, the only 
difference between the center and the layer of the jet is 
the bulk Lorentz factor which is determined for the spine 
(r spine) by modeling the BL Lac SED, while for the layer 
(Y layer) by requiring that the debeamed BL Lac match the 
FR I distributions in the L ext — L Q and L r — L Q planes. 
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Fig. 8. Debeaming trails in the radio-optical luminosity 
plane, for the case of the two-velocity jet. The five curves 
correspond to average SED (1,2 and 3), Mkn 421 (M) and 
PKS 0735+178 (P). For the average SEDs T spme = 15 
and ri a y er — 2; for Mkn 421 r spirie = 20 and Ti ayer = 1.2; 
for PKS 0735+178 T spme = 16 and T layer = 1.5. Circles 
correspond to the predicted optical-radio luminosity for 
different angles of sight (top to bottom: d=l/T, 10°, 30° 
and 60°). 
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Fig. 9. Debeaming trails for the optical emission of 
Mkn 421 and PKS 0735+178 in the two-velocity jet sce- 
nario. The circles connected by the vertical lines corre- 
spond to the sources observed at angles of (top to bottom) 
1/r, 10°, 30°, 60° (filled black circles) and 90°. The values 
of Ti ayer (1.2 for and 1.5 for Mkn 421 and PKS 0735+178, 
respectively) are chosen in order for the luminosity at 60° 
to correspond to the median value for each bin of extended 
radio power. 



The monochromatic intensity emitted by the jet is 
therefore calculated as 

IW) = C(^'WW) + SlayeriWiW Slayer), 

where /' is the comoving intensity. 

The predicted luminosity trails for the two specific BL 
Lacs are shown in Fig. |^. Values of Ti ayer are set to 1.2 and 
1.5 for Mkn 421 and PKS 0735+178 respectively, so that 
the point of each trail corresponding to the angle 9 — 60° 
falls approximatively onto the median of the FR I optical 
core luminosity in each L ex t bin (Fig. |^). 

By using the same value of Ti ay er this two-velocity 
model satisfactorily predicts the properties of the de- 
beamed counterparts of LBL and in particular it repro- 
duces the FR I location in the L r — L Q plane. 

Conversely this picture can not account for the ob- 
served optical-radio properties of debeamed HBL. The 
luminosity of objects seen at 9 = 60° is close to the lower 
limit of the FR I region in the optical, but still one order 
of magnitude fainter in the radio. We must stress how- 
ever that the extended radio powers of HBLs correspond 
more closely to the range covered by the B2 radio galaxies. 
Clearly, any firm statement on this issue must await for 



the analysis of the nuclear properties of the B2 sample, 
but the extrapolation of the 3CR radio-optical correlation 
does not match the debeamed predicted luminosities of 
HBL. This result does not depend on the specific value of 
Flayer adopted since, as already discussed in Sect. [|, the 
HBL trails run almost parallel to the radio-galaxies locus. 

A further modification of this model is thus required 
for the HBL unification. Without altering the comoving 
spectra of the two components, the simplest change is 
to assume a lower Doppler factor for spine in the radio 
emitting region, as might be the case if the flow slows 
down between the optical and the radio emitting sites. 
This would increase the initial slope of the debeaming trail 
which would rapidly reach the FR I region. 

5.1. Constraints from the X-ray observations 

The limited angular resolution makes the analysis of X-ray 
observations of FR I sources less straightforward than in 
the radio and optical bands. In particular it is necessary 
to disentangle any non thermal nuclear radiation from the 
often dominant emission of the hot gas associated with the 
galactic corona and/or galaxy cluster. Nonetheless they 



8 



Chiaberge et al.: FR I - BL Lacs unification 



N 



o 

l-J 




22 24 26 28 
Log L x [erg s _1 Hz -1 ] 



Fig. 10. Debeaming trails in the optical- X-ray luminos- 
ity plane for the average BL Lacs SED. Two-velocity jet 
with r 



spine 



15 and Ti ayer — 2 for all SED. Circles are 
FR I, squares and triangles represent X-ray and radio se- 
lected BL Lacs, respectively. Grey scales correspond to 
bins of extended radio luminosity, as for Fig. |l|. Open tri- 
angles and squares are sources without measurement of 
their extended radio power. Open circles correspond to 
the predicted optical- X-ray luminosity for different angles 
of sight (top to bottom: 9=1/T, 10°, 30° and 60°). 



provide useful constraints to high energy nuclear emission 
of these radio-galaxies. 

We can test the validity of the two-velocity jet sce- 
nario by considering also this X-ray emission. In Fig. [lO] 
we report the debeaming trails in the Lx — L Q plane for 
the average SEDs assuming the same jet parameters as 
before: the predicted powers in both bands appear to be 
consistent with the observed properties of radio galaxies, 
supporting the presence of a less beamed plasma compo- 
nent (layer) dominating the emission in the parent popu- 
lation also in the X-ray band. Conversely, there is no need 
for a different amount of beaming in these two bands. This 
is somehow reassuring, as optical and X-ray are believed 
to originate co-spatially. 

In addition, X-ray data can be used to define the loca- 
tion of radio-galaxies also in the broad band spectral in- 
dices plane. The spectral characteristics of blazars are of- 
ten represented in the plane defined by a ro (5GHz-5500A) 
and a oa; (5500A-l keV). It is therefore worthwhile to de- 
termine how relativistic beaming affects the position of 
the objects also in this plane. While an approximated re- 
lation between the BL Lacs and F R I broad band spectral 
slopes is derived in Appendix A: in the case of constant 



local spectral indices, changes in the local spectral slopes 
are properly taken into account. 

In this plane (Fig. [Tl] ), as already well established, 
HBL and LBL occupy the left (i.e. flatter a ro ) and the 
top-center (i.e. steeper a ox ) regions, respectively and their 
different position is acco unted for by their different SEDs 
(e.g. Fossati et al. 199? ), i.e. reflects the position of the 
peak of the emission. The FR I region is instead well de- 
fined at the center of the diagram. The debeaming trails 
for Mkn 421 and PKS 0735+178 are also shown f| the 
empty circles correspond to 9 — 1/T, 10° 30°, 60° in the 
case of the two components model, while the two aster- 
isks represent each source as observed at 60° in the case 
of a single emitting region. PKS 0735+178 falls in the ra- 
dio galaxy region for Ti ayer = 3 or less, while Mkn 421 
does not intersect this area either in the single or in two 
component models, confirming the results of the analysis 
presented above. 



6. Summary and conclusions 

With the aim of exploring the viability of the unification 
scenario between (HBL, LBL) BL Lacs and FR I radio 
galaxies we have compared their nuclear emission in the 
radio, optical and X-ray bands. 

We have firstly considered these spectral regions sepa- 
rately, comparing the nuclear emission of the two classes 
of objects for similar extended radio power. As the core 
radiation of BL Lacs is enhanced by relativistic beaming, 
we derived the bulk Lorentz factors requested to account 
for the observed distribution. The values of T thus in- 
ferred are not compatible with the higher bulk velocities 
requested by theoretical arguments, such as the pair pro- 
duction opacity and the spectral modeling of the SED of 
BL Lacs. 

We then examined the core emission of three samples 
in the L r — L a plane. In the frame of the simplest one-zone 
emission model, we calculated debeaming trails of the BL 
Lac broad band emission as predicted by the relativistic 
transformation for an increasing angle of sight. We found 
that the model does not account for the observed spectral 
properties of FR I, as expected from the above inconsis- 
tency of the Lorentz factors. 

The simplest and rather plausible hypothesis to ac- 
count for this discrepancy within the unification scenario 
is to assume a structure in the jet velocity field, in which a 



5 Note that again the two trails differ because of the differ- 
ent SEDs: for Mkn 421 the peak of the synchrotron component 
is between the optical and the X-ray bands and therefore the 
effect of debeaming is, initially (i.e. for increasing 9) to steepen 
both a ro and a ox ; when the Compton peak enters the X-ray 
band, a ox flattens. Instead, for PKS 0735+178 the synchrotron 
peaks between the radio and the optical, and the 1 keV flux is 
due to Compton and, for the range of angles of sight consid- 
ered here, the change in the spectral indices turns out to be 
monotonous. 
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Fig. 11. Debeaming trails in the a ro — a ox plane for Mkn 
421 and PKS 0735+178. The jet parameters are the same 
as in Fig. |[ The black filled symbols correspond to the 
observed BL Lac (9 = 1/F), while the empty circles rep- 
resent the predicted position for different angles of sight 
{9 = 10°, 30° and 60°) in the frame of a two- velocity 
jet model. The two asterisks indicate sources observed at 
9 = 60° in the case of a single emitting component. 

fast spine is surrounded by a slow layer. Note however that 
the slower jet component must be relativistic in order to 
explain the anisotropic radiation of radio galaxy cores (e.g. 
Capetti & Celotti 1999). The observed flux is dominated 



by the emission from either the spine or the slower layer, 
in the case of aligned and misaligned objects, respectively. 

Interestingly, the existence of velocity structures in the 
jet has been suggested by various authors (Komissarov 
199C, Laing 1993| ) in order to explain some observed prop- 



erties of FR I (and FR II) jets, such as the structure of the 
magnetic field in FR I which appears to be longitudinal 
close to the jet axis and transverse at the edges. Swain 
et al. |1998| obtained VLA images of 3C 353 (an FR II 
with straight jets), finding that a model consisting in a 
fast relativistic spine (/3 > 0.8) plus a slower outer layer 
{(3 < 0.5, but still relativistic in order to produce the ob- 
served jet-counterjet intensity asymmetry) could account 
for the apparently lower emissivity near the jet axis. Simi- 
lar behaviours have been inferred for the two low luminos- 



ity radio galaxies M 87 (Owen et al. |1989| ) and B2 11 44+35 
(Giovannini et al. 1999] ). Furthermore, Laing et al. (1999) 
showed that the jet asymmetries in FR I can be explained 
by means of a two-speed model. As a consequence, they ar- 
gued that the lower velocity component dominates in the 
cores of the edge-on sources, while the fast spine emission 



dominates the end-on ones. This possibility might be also 
supported by recent numerical simulations of relativistic 



jets (Aloy et al. |2000[ ). 

The same indication has been f ound through differ- 
ent approaches. Capetti & Celotti 199S| reveal a trend 
in the radio galaxy /BL Lac relative powers with the line 
of sight, which is consistent with a slower (less beamed) 
component dominating at the largest angles. Capetti et 
al. 2000 consider the same issue by examining the more 
detailed SED of five radio galaxies and consider their 
beamed counterparts. They found that while the spec- 
tral shapes of 3C 264 and 3C 270 can be reconducted to 
those of BL Lacs, the required ratio of beaming factors, 
i.e. ^BLLac/^FR i ~ 10 — 100, implies that the correspond- 
ing BL Lacs would be ovcrluminous. The inclusion of a 
slower (less beamed) jet component seems to be a plausi- 
ble explanation. 

We found that Lorentz factors of the layer ri aver ~ 2 
can account for the unification of FR I (of the 3CR) with 
LBL and intermediate luminosity BL Lacs. Instead the de- 
beaming trails for the lowest luminosity HBL do not cross 
the FR I region in the L r — L a plane. While the HBL 
behavior should be compared with that of radio galaxies 
with which they share the extended radio power (e.g. those 
of the B2 catalogue), our simple two-component jet model 
could not account for the observed properties if the cores 
of such low-power FR I radio galaxies lied on the extrapo- 
lation of the 3CR radio-optical correlation. The properties 
of such weak sources can be instead reproduced if their ra- 
dio emitting region is less beamed than the optical one, 
as could be expected if the jet decelerates after the higher 
energy emitting zone. 

Finally, the presence of velocity structures in jets of 
course affects the number counts of beamed and unbeamed 
sources: for example, the lack of BL Lacs in clusters (Owen 
et al. 1996] ) could be attributed to values of typical bulk 
Lorentz factors higher than those derived from statistical 
arguments (Urry et al. 1991). Intriguingly, the very latter 
authors had to require a wide distribution of Lorentz fac- 
tors to account for the number densities of FR I and BL 
Lacs in the radio band. 

Much has still to be understood on the dynamics 
and emitting properties of relativistic jets. Multifrequency 
studies of the nuclear properties of beamed sources and 
their parent populations and their comparison - accord- 
ing to unification scenarios which are well supported by 
other independent indications - constitute a new and pow- 
erful tool to achieve that, both for well studied individual 
sources as well as complete samples. Near IR observations 
by HST, mm data and higher resolution and sensitivity 
by Chandra in X-rays will further open this possibility. 

Concluding, the radio, optical and X-ray nuclear emis- 
sion of FR I and BL Lacs strongly indicate the presence 
of a velocity structure in the jet if indeed these sources 
are intrinsically identical. In other words, by considering 
the indications of trends in the SED of blazars emerged 
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in the last few years (Giommi & Padovani 1994, Fossati 
et al. 1998) together with the constraints derived from 



their unification with radio galaxies, it appears that the 
phenomenology of these sources is characterized and de- 
termined by differences both in the intrinsic SED it and in 
beaming properties. 
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Appendix A: Debeaming and the broad band 
spectral slopes 

In the case of a one-component model, and under the as- 
sumption that the local spectral indices are constant, we 
can derive the transformation law for the change of the 
spectral slope due to relativistic beaming. If the flux den- 
sity in the frame comoving with the emitting region is F' v , , 
the observed one is 

where <5 bjcct is the beaming factor of the same object 
for different lines of sight (i.e. observed as BL Lac or as 
radio galaxy). Substituting these transformations in the 
definition of the broad-band spectral index a\2 (where 1 
and 2 refer to radio, optical or X-ray flux), one obtains 
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which can be written as 
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where ai and a 2 are the local spectral indices. 



